Nonlinear bias dependence of spin-transfer torque from atomic first principles 
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We report first-principles analysis on the bias dependence of spin-transfer torque (STT) in 
Fe/MgO/Fe magnetic tunnel junctions. The in-plane STT changes from linear to nonlinear de- 
pendence as the bias voltage is increased from zero. The angle dependence of STT is symmetric at 
low bias but asymmetric at high bias. The nonlinear behavior is marked by a threshold point in 
the STT versus bias curve. The high-bias nonhnear STT is found to be controlled by a resonant 
transmission channel in the anti-parallel configuration of the magnetic moments. Disorder scattering 
due to oxygen vacancies in MgO significantly changes the STT threshold bias. 

PACS numbers: 72.25.Ba, 72.10.Bg, 85.75.-d 



I. INTRODUCTION 

Magnetic tunnel junction (MTJ) has attracted great 
attention due to its importance in magnetic random- 
access memory (MRAM) and read-sensor technology. A 
basic MTJ is made of two ferromagnetic (FM) layers 
sandwiching a thin insulating material. Digital informa- 
tion is coded by magnetic moments of the FM layer being 
in parallel or anti-parallel configurations (PC or AFC). In 
commercial MRAMs, switching between PC and APC is 
achieved by external magnetic fields. An emerging trend 
is to switch by spin-transfer torque (STT) predicted 
theoreticall}^ and demonstrated experimentally^Ti^ where 
a spin-polarized current transfers angular momentum to 
a FM layer causing magnetic moment reversal. STT 
holds great promise to simplify the MRAM structure, 
down-scale device size, and reduce power consumption. 

One of the most important issues concerning STT is 
its dependence on external bias voltage VJ, that drives 
the spin-polarized current in the first place. Understand- 
ing this issue is rather difficult not only because torque 
is a vector—"— but also because STT operates at non- 
equilibrium due to the flow of spin currents. A physi- 
cal picture concerning the bias dependence of STT, from 
the linear to nonlinear bias regime, is yet to be estab- 
lished from atomic first principles. Theoretically, Vf, 
dependence of STT was investigated by model analy- 
sis such as the tight-binding mode l^°i^^ and free-electron 
modelii^rii More recently, Heiliger and Stiles^^ analyzed 
STT for an MgO based MTJ based on equilibrium elec- 
tronic structure obtained from density-functional the- 
ory (DFT). The in-plane STT was calculatedi^ up to 
Vf, = ±0.5 V and showed essentially a linear depen- 
dence on Vb, and the out-of-plane STT was found to be 
quadratic in 14, which is consistent with tight-bindingiS, 
and experimental resultsi^ii^ 

Existing calculations provided valuable understanding 
of STT at the small Vf, regime. To establish a complete 
picture, it is important to investigate STT at higher bias 
from atomic first principles. Indeed, recent experimental 
datai^ already showed a strong nonlinear Vf, dependence 



of in-plane STT, and observable non-quadratic out-of- 
plane STT beyond Vf, w 0.2V. It is the purpose of this 
work to report a first-principles analysis of STT in the 
most technologically important Fe/MgO/Fe MTJ from 
low- to high-bias regimes, focusing on its nonlinear Vf, 
dependence. We found that the in-plane STT changes 
from linear to nonlinear dependence as Vf, is increased 
from zero. The angle dependence of STT is symmetric 
at low bias but asymmetric at high bias. The nonlin- 
ear behavior is marked by a threshold point in the STT 
versus bias curve which is controlled by a resonant trans- 
mission channel in the APC current. We also found that 
the nonlinear STT is affected by atomic defects— oxygen 
vacancies in MgO— hence, in principle, can be tuned by 
interface engineering. Results are compared to the cor- 
responding experimental data. 

This paper is organized as follows. In Sec. II, we give 
the details of our calculation based on Keldysh nonequi- 
librium Green's function (NEGF) and wave-function- 
matching function method. In Sec. Ill, we present our 
results on Fe/MgO/Fe(001) MTJs with clean and disor- 
dered interfaces. Section IV is our summary. 



II. ELECTRONIC STRUCTURE AND 
TRANSPORT CALCULATION 

The MTJ we consider consists of an MgO barrier sand- 
wiched by two semi-infinite Fe leads shown in FigUJ The 
system is periodic in the x-z plane and current flows 
along y corresponding to the (001) material-growth di- 
rection. A very small lattice mismatch between Fe and 
MgO is neglected by fixing the interfacial atoms at their 
bulk bcc positions. Our calculation is from first principles 
where DFT is carried out within the Keldysh nonequilib- 
rium Green's function (NEGF) method^^ Here, DFT 
determines the Hamiltonian of the open device struc- 
ture, and NEGF determines the nonequilibrium statis- 
tics of the device operation and the nonequilibrium den- 
sity matrix. The NEGF-DFT is solved self-consistently 
under finite bias during current flow. For devices hav- 
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FIG. 1: (Color online) Sketch of a Fe/MgO/Fe(001) MTJ 
having seven MgO layers. The magnetization of the left Fe 
lead (Ml) is fixed and that of the right Fe lead (M2) is free 
to rotate. M2 is pointing along the z axis; Mi lies in the x-z 
plane, the angle Q is the angle between Mi and M2. The red 
and blue grids in the scattering region (MgO) denote the O 
and Mg atoms, respectively. There may be some interfacial 
disorder (oxygen vacancy, cyan grid) at the Fe/MgO inter- 
faces. The applied bias eV(,=/iii — [il, where [ir and are 
chemical potentials of the right and left leads. 



ing oxygen vacancies (OV, the most energetic favorable 
defects), we further apply the nonequilibrium- vertex- 
correction (NVC) theorjii^ for disorder averaging at both 
the nonequilibrium-density-niatrix level and transport- 
calculation level. The NEGF-DFT-NVC formalism al- 
lows one to self-consistently calculate quantum trans- 
port properties from atomic first principles without phc- 
nomcnological parameters. Wc refer interested readers 
to the original workJ^ 

Moreover, The self-consistent nonequilibrium poten- 
tials were used as input to a TB-MTO wave-function- 
matching calculation. The scattering wave functions of 
the whole system were obtained explicitly. Because wc 
cannot embody the nonequilibrium NVC into the out- 
of-planc STT calculations at present, a supercell is used 
to assess out-of-plane STT in the presence of interfacial 
disorder. In the calculations, an 800 x 800 A;-mesh is 
used to sample the two-dimensional Brillouin zone (BZ) 
in the x-z plane (see Fig. [ij to ensure accurate conver- 
gence. Other numerical details are similar to those of 
Ref. [HandHS 

To check our electronic structure of Fe/MgO/Fe, we 
calculate the zero-bias TMR for 7 monolayers (ML) MgO 
barrier, which is 5100% with the perfect Fe/MgO inter- 
face while decreases dramatically to 350% (78%) with 
3% (7%) OV at both interfaces. The TMR of a clean 
junction is consistent with the published first-principles 
calculation;^ while the TMR of junctions with vacan- 
cies are in the range of experimental measurement 1^^— 
The results indicate that wc have got the right electronic 
structure. 



III. BIAS DEPENDENCE OF STT 

Before presenting results, let's consider some general 
trends. Neglecting any spin dependent scattering in the 
tunnel barrier, the in-plane STT can be obtained from 
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FIG. 2: (Color online) (a) STT of ideal Fe/MgO(7L)/Fe with 
angle 6 = 90° between Mi and M2. Black squares and red 
circles are out-of-plane and in-plane STT calculated by the 
wave-function-matching method, respectively, the blue up- 
triangle is the in-plane STT calculated by Eq. (1). Inset: spin 
currents in PC and APC. (b) and (c): angular dependence of 
STT at Vt = -0.4, and -0.75 V, respectively, (d) and (e): 
angular dependence of TMR at Vt — —0.4, and— 0.75 V, re- 
spectively; the red dash lines are standard cosine-dependence 
TMR according to Ref. ^ 



conservation of spin current^ 
= 



M2 X (Ml X M2) (1) 



where 6 is the relative angle between magnetization Mi 
and M2 (see Fig. H]). li'\o) = ^[I^{0) -1^(0)] is 

the spin-polarized current density in PC and /i''''(7r) = 
^[/^(tt) - I^{Tr)] is that of APC. By structural symme- 
try, /i*'' (0) should be an odd function and (tt) an even 
function of Vb- Hence, at not too large Vb, we can write 
/i^'(O) « 7iVfc and /i'*''(7r) « 72 V^^, where 71 and 72 are 
constants. From these considerations, Eq. ^ suggests 
Iz^\o) to dominate at small Vb (positive Vb) while 
/i^''(7r) to dominate at higher Vb- Hence, Tm should start 
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from a linear dependence and turn into a nonlinear de- 
pendence as Vt is increased; the turning point defines a 
"threshold" bias whose value is determined roughly by 
the material-specific constants 71,2- At low bias, the out- 
of-plane STT demonstrates a simple quadratic relatiorti^ 
that we can write as Tj_ = j±_V^. 

After the NEGF-DFT self-consistent calculation is 
convergedr^ spin currents and can be calculated 
by NEGFi^ and the in-plane STT T|| obtained from Eq. 
(P. The out-of-plane STT T± (also T||) can be obtained 
from a scattering-wave-function approach following Rcf. 
[19I . Figure [2] shows plots of the calculated T|| and T± ver- 
sus Vb for a perfect Fe/MgO/Fe MTJ (no disorder) hav- 
ing seven layers of MgO, where magnetization Mi _L M2, 
i.e., 6 = 90° (see Fig. [T|). Seven-layer MgO is used be- 
cause it is the thickness of experimental devicei^ with 
which we shall compare results (see below). The blue 
(up-trianglc, T||) and black (solid square, T±) curves were 
calculated from the scattering wave functionsJ^ The red 
curve (solid circle, T| | ) was calculated from Eq. ([T]) , which 
agrees with the blue curve very well at least for lower 
At lower bias up to ±0.5V, Tn is linear in Vb and it 
is in agreement with Reffisl As Vb increases, Tn becomes 
gradually nonlinear and a clear threshold at Vb ~ 0.65V 
(e.g., turning point of the curves) is seen, as expected 
from the general considerations discussed above. 

The inset of Fig. [H^a) plots the spin-current density 

versus Vb- Spin-polarized current density /i*^(0) is in- 
deed linear in Vb whose slope is 71 « 35/iJV^^m^^. Spin- 
polarized current density /i''^ (tt) is, as expected, an even 
function of Vb and, by fitting to a quadratic form 72 Vj,'^, 
we found that 72 w S/^JV-^m-^ for Vb < 0.5 V. The 
small value of 72 makes li^^ (tt) to have a very weak Vb 
dependence, until its sudden increase at Vb ^ 0.6 V, re- 
sulting in the threshold of Tn at ^0.65 V. 

The out-of-plane STT shows quadratic behavior at the 
lower bias with j± « — 13//JV^^m^^, and it is very close 
to the published calculated value of — 14/iJV^^m^^ with 
a SIX- ML MgO barrier Quantitatively, Tj_ < T\\ for 
Vb < 0.5 V but reaches a similar scale when Vb > 0.5 V. 

Not only the bias dependence of STT [see Fig. [H^a)] is 
changed by the sudden increase of APC current /i'*-'(7r), 
the angular dependence of STT is also modified. As 
shown in Fig. [^b), the angular dependence behaves as 
sin^ at a small bias Vb = —0.4 V. It however deviates 
from the sine dependence at a higher bias Vb = —0.75 
V [sec Fig. IH^c)]. Experimental evidence of asymmetric 
angular-dependent STT at higher bias had been already 
observed before^i^ 

Can the angular dependence of in-plane STT be under- 
stood by the angular dependence of the charge current 
or TMR? We give the angular dependence of TMR at a 
small bias Vb = —0.4 V [see Fig. [D^d)], and a large bias 
Vb = —0.75 V [see Figl2je)]. The angular dependence of 
TMR at low bias follows the cosine function^^ and devi- 
ates from the cosine function at high bias. At high bias, 
the in-plane STT shows an asymmetry and the peak is 
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FIG. 3: (Color online) Energy dependent spin polarized trans- 
mission coefficient of ideal Fe/MgO(7L)/Fe MTJ at 14 = 
±0.75V. Insert (a); the k|| resolved transmission hot spot 
of the t -channel in APC at energy E^EfH- 0.36eV. Insert 
(b): the energy dependent transmission of a bright point A 
on the ring of the insert (a) at ka:=kj,=0.278A~^ . Insert (c): 
angular dependent in-plane STT at resonant point A kx=ky= 
0.278A"\ and normal point B k^^ky= -1.624A"^ of the in- 
sert (a); therein. To = ^kQ^^m^^ sets the unili^^ 



close to the parallel side, while the angular-resolved TMR 
decreases always faster than cosine function as the mag- 
netization goes from parallel to antiparallel. Notice here, 
we have negative TMR at high bias. So, the angular de- 
pendence of TMR cannot explain our calculated angular 
dependence of in-plane STT. 

It is the sudden increase of APC current /i'*''(7r) at a 
higher bias Vb = 0.6 V [inset of Fig. EJa)] that led to the 
nonlinear bias dependence of STT. To understand the 
microscopic origin of this sudden increase, we plot spin- 
resolved transmission coefficient versus electron energy E 
(relative to the Fermi energy Ep) in Fig. [3l calculated at 
Vb = 0.75V. The t or J, denotes the majority or minority 
spin channel in the left Fe lead. For APC, the t channel 
of the left lead transmits to the J, channel of the right 
lead through the MgO. Figure (3] shows that there is an 
abrupt increase in f channel transmission at about E = 
0.3 eV. For a symmetric MTJ, roughly half Vb is dropped 
at each Fe/MgO interface. Hence the abrupt increase of 
transmission at E w 0.3 eV should contribute to APC 
current when the bias Vb — 0.6 V. This is indeed what is 
found [inset of Fig Elja)]. We conclude that the sudden 

increase of li"^ (tt) is due to the abrupt increase of the t 
channel transmission, which, importantly, is also due to 
bias-dependent transport features as we explain now. 

The inset (a) of Figl3] plots the transmission hot spot 
of APC(t) in the two-dimentional (2D) BZ, which is the 
transverse momentum (k||) resolved transmission coef- 
ficient. There is a clear bright ring surrounding the F 
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FIG. 4: (Color online) STT of Fe/MgO(7L)/Fe with oxy- 
gen vacancy disorder versus bias. LxRy denotes that there 
are x and y OVs at the left and right interface, respectively. 
N0.03O0.97 denotes that 3% oxygen atoms at both interfaces 
are replaced by nitrogen atoms. 



point. The abrupt increase of the APC t channel trans- 
mission is due to this bright ring. Focusing on one partic- 
ular k point on this ring and scanning the electron energy, 
a single resonant peak is discovered in the transmission 
as shown in the inset (b) of Fig. [31 

Importantly, this sharp resonance is different from the 
well-known resonance states on Fe surfaces.— In partic- 
ular, we found that this sharp APC resonance is only 
established at high bias and is located between the Ai 
band in the left Fe and the A5, A2/ bands in the right Fe. 
The contribution of A5 and A2/ bands to the total trans- 
mission is k|| dependent, giving rise to the ring in the 
hot spot. A higher bias significantly changes the shape 
of the tunnel barrier: we checked that by artificially fiat- 
tening the barrier potential from the biased slop shape, 
this sharp APC resonant channel disappears. 

Actually, the existence of resonance is not only respon- 
sible for the sudden increase of STT but also responsible 
for the deviation of angular-dependent STT as shown in 
the inset (c) of Fig. [3] When at a resonant ring such as 
A point {kx = ky = 0.278 A~^), the angular- dependent 
in-plane STT deviates from a sine function and follows 
similar behavior as that in Fig. [H^c). When far from 
the resonant ring such as B point (fcj. = ky = —1.624 
A"^), the angular dependent in-plane STT follows a sine 
function as well. Under most situations, at high bias, 
the resonance dominates the in-plane STT. Without res- 
onant states, the Bloch wave should decay very fast in 
the tunnel barrier, the multiple scattering should be neg- 
ligible, and Eq. [1] should hold in this case. However, 
in the presence of resonant states, the wave function of 
these states will decay little in the barrier and Eq. [T] is 
not valid anymore. 

Having understood the microscopic physics behind the 
nonlinear bias dependence of STT, we note that the ob- 



tained threshold bias at H ~ 0.6 V for the ideal junction 
is higher than that reported experimentally^^ which was 
about 0.3 V. Since experimental devices are never ideal, 
an investigation of defect scattering and its effect on STT 
is warranted. There are several possible defects such as 
oxygen vacancy in MgO, interfacial roughness^^ material 
imperfections, etc. It was reported experimentally that 
due to compressive strain during crystal growth, oxy- 
gen vacancies (OV) are accommodated in MgOi^ Both 
theory2£ and experiment!^ have shown that existence of 
OV can drastically infiuence spin-polarized transport and 
its effect on STT is therefore expected. 

We calculated the in-plane STT from Eq. (P) by 
putting OVs in the interfacial MgO layer immediately 
adjacent to Fe. An alloy model Oi-^i-Va^; is used, where 
Va stands for vacancy and x is the OV concentration. 
The disorder average is carried out by the NVC theoryi^ 
mentioned above. Results are plotted in Fig. 2] versus 
Vf); and several observations are in order. 

First, at low VJ,, T\\ is linear for both disordered {x ^ 0) 
and clean samples {x = 0) — hence, the linear regime is 
not qualitatively altered by OVs. Second, both sym- 
metric junctions where the same OV concentration is on 
both Fe/MgO interfaces (black curve with solid squares), 
and asymmetric junctions where OVs only exist at one 
Fc/MgO interface (red curve with solid circles), yield al- 
most the same T\ | . In particular, its slope r| | = dT\ | /dV ~ 
—5.5 X IO^'^Tq. This value is only slightly larger than 
that of the ideal junction, which has t\\ = —4.7 x 10^'^ Tq. 
Hence, a small amount of OV (e.g. 3%) does not alter 
T|| significantly. These results should be compared to 
the experimental value of the slope r|| = —3.2 x 10^"^ Tq 
measured at Vb ~ ±0.3 VJ^ Giving the uncertainties 
in comparing our atomic structure to the experimental 
one, the quantitative consistency of our results is quite 
satisfactory. Third, in Ref. [2^ it was shown that fill- 
ing the OVs with nitrogen atoms can significantly reduce 
disorder-induced diffusive scattering and drastically in- 
crease the TMR ratio. Here, filling the OVs with N using 
the alloy model Oi-ajN^;, we found that the Tn is some- 
what decreased, namely, its slope r|| = —1.9 x IO^'^Tq 
(blue curve with up-triangles). 

Similar to the ideal junction, systems with OV also 
show abrupt increase of T|| at higher bias, which, we 
found, is also due to the quantum resonance in the 2D BZ 
discussed above. Different from the ideal case, existence 
of OVs shifts the threshold bias of T| | to lower values thus 
bringing it closer to the experimental result 1^ As shown 
in the inset of Fig. |4l the effect of OV on | can also be 
understood by investigating transmission of the t channel 
in APC. 

The OVs induce two new effects. First, interfacial OVs 
broaden the area having high transmission values in the 
hot spots and also suppress the transmission peak from 
1.0 e'^/h to about 0.35 e^/h (at 3% OV on both inter- 
faces). The overall effect is to enhance T|| slightly. Sec- 
ond, interfacial OVs shift the resonant peak in the t chan- 
nel of APC to lower energy and, as a result, the threshold 
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bias is now shifted to lower values. 

Remarkably, from Fig. 2] one can determine which 
Fe/MgO interface dominates STT. When OVs exist at 
both interfaces (black curve with solid squares), the 
threshold bias shifts to lower values from those of the 
clean sample for both Vt polarity. When OVs exist only 
at the right Fe/MgO interface (red curve with solid cir- 
cles), the threshold bias shifts to lower value only in the 
negative VJ,. When OVs are filled by N atoms (blue curve 
with up-triangles) , the threshold bias shifts to higher Vj, 
for both polarities. 

Quantitatively, the threshold bias is ^0.45V in the 
presence of 3% OV at interface. It is consistent with the 
experimental value '-^0.4 V (see Ref.— ) with a 1.25-nm 
(6-ML) MgO barrier. A more recent experiment mea- 
sured a threshold bias '^0.2 V (see Ref4) with a 1.0-nm 
(5-ML) MgO barrier. The experimentally measured re- 
sistance is 1.5 f2/im^ also close to our calculated value for 
a junction with a 5-ML MgO barrier. 

The effects of OVs on the out-of-plane STT are also 
very important. We have examined a junction using the 
wave-function-matching method. The calculation was 
done in a 7 X 7 lateral supercell containing five ran- 
domly distributed OVs at each Fc-MgO interface (thus 
X - 10%). For Vb = 0.5 V and £■ = Sf + 0.25 eV, the cal- 
culated in-plane and out-of-plane STT are 20.30 x IO^'^Tq 
and 7.32 x IO^'^Tq, respectively. These values are larger 
than that of the ideal junction which are 1.06 x IO^^Tq 
and 0.94 x IO^^Tq for in-plane and out-of-plane STT, re- 
spectively. The larger in-plane and out-of-plane STT in 
the 10% OV sample is attributed to the resonant tun- 



neling discussed above. These results indicate that the 
intcrfacial OV disorder is not detrimental to the value 
of the out-of-plane STT. In comparison, interface rough- 
ness disorder quenches the out-of-plane STT as reported 
in Ref. HI. 



IV. SUMMARY 

In summary, the bias dependence of STT was inves- 
tigated from atomic first principles, covering both low- 
and high-bias regimes. STT shows weak and linear bias 
dependence at small Vt but strong and nonlinear depen- 
dence at higher Vt- The angle dependence of STT is 
symmetric at low bias but asymmetric at high bias. The 
nonlinear behavior is marked by a threshold bias in the 
STT versus bias curve and is controlled by a resonant 
transmission channel in AFC. Very importantly, disorder 
scattering by small amount of oxygen vacancies in MgO 
lowers the STT threshold bias but not the STT value, 
suggesting that the nonlinear STT can be tuned by clever 
interfacial engineering as also reported in Ref. |30| where 
crystal-growth temperature changes the OV contents. 
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